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Abstract 

The fragmentation of the ground state of a repulsive condensate immersed into a double-trap 
potential is found to be a general and critical phenomenon. It takes place for a given number 
of bosons if their scattering length is larger than some critical value or for a given value of the 
scattering length if the number of bosons is above some critical number. We demonstrate that the 
geometry of the inner trap determines these critical parameters while the number of the fragments 
and the fraction of bosons in the various fragments can be manipulated by the outer trap. There 
is also a maximal number of bosons for which the ground state is fragmented. If this number is 
exceeded, the fragmented state becomes a very low-lying excited state of the condensate. This 
maximal number of bosons can be substantially manipulated by varying the inner and outer traps. 
To study three-fold fragmentation we have chosen a potential well with two barriers as the inner 
trap and embedded by two types of outer ones. A many-fold fragmentation is also addressed. 

PACS numbers: 03.75.Hh,03.65.Ge,03.75.Nt 
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I. INTRODUCTION 



The general mathematical formulation of the condensation phenomenon for an ideal gas 
in equilibrium has been given by Penrose and Onsager Q in 1956. They considered a sys- 
tem of N interacting bosons and its reduced one-particle density matrix. The eigenvalues 
of this matrix are called occupation numbers and the eigenvectors are referred to as natu- 
ral orbitals. If an ideal gas of free bosons forms a BEC, then only a single natural orbital 
is macroscopically occupied. The extension of this concept on to several macroscopically 
occupied orbitals is the basis of the fragmentation phenomenon. The BEC is called frag- 
mented if several natural orbitals have macroscopic occupation numbers Q. Although these 
original definitions have been formulated in the thermodynamic limit N —>■ oo, the reduced 
one-particle density matrix can also be applied to study condensation [3] and fragmentation 
^| in finite-N bosonic systems. 

The first measurements of an interference between two expanding and overlapping con- 
densates j^| have stimulated a great interest in the theoretical studies of fragmentation. An 
initial quantum state of such a system is supposed to be two-fold fragmented i.e. two spa- 
tially separated orbitals are macroscopically occupied. The theoretical studies on relative 
phases p, bllg between the fragments and related questions on dynamical evolutions of 
this state 9l llO | have been a subject of numerous discussions. A natural extension of these 
effects to an array of multiple wells with a many-fold fragmented initial state has been 
initiated by the recent experiments in optical [l^ and hybrid traps jl^ . 

In contrast, the fragmentation of the ground state of BEC has been predicted only for 
a few systems where it is enforced by the spatial or spin symmetry of the systems. In 
particular, the ground state of an attractive condensate in a perfect ring is found to be 
fragmented [14]. The rotational symmetry of the corresponding Hamiltonian permits as 
natural orbitals only plane waves 15]. The many-body wave-function of this system starts 
to differ from a function composed by plane-waves for any non-zero attractive interaction 
and the respective one-particle density matrix having plane-waves as natural orbitals must 
exhibit several non-zero eigenvalues. This eventually leads to fragmentation and its origin is 
spatial (rotational) symmetry. Another example is an attractive BEC in a symmetric double- 
well potential where the ideal symmetry of the potential 16 1 causes the fragmentation of the 
ground state. The stability of all these fragmented states with respect to a small asymmetric 
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perturbation is still an open question. There are indications that these symmetry related 
fragmentations disappear upon symmetry breaking. Fragmentation can also take place if 
bosons have additional internal degrees of freedom. An example of such a system is a 
system made of Bose particles with non-zero spin, a so-called spinor condensate. In spinor 
condensates different spin components may have different spatial extensions in the presence 
of an external magnetic field. Indeed, the ground state of Spin-1 Bose gas in a uniform 
magnetic field was proved to be fragmented jl?! ]. 

The single condensate picture and its mean-field description via the Gross-Pitaevskii (GP) 
equation has been a very successful approximation and can explain many experiments, see, 



e.g., Refs. 



la . Il9| and references therein. However, this mean-field is incapable by definition 



to describe fragmentation, since only one orbital is involved. Recently, a more flexible mean- 
field approach allowin g fo r bosons to reside in different orthonormal one-particle functions 
has been formulated [16]. This intrinsic ability to describe fragmentation makes it very 
attractive for theoretical investigations and predictions. In the framework of this best mean- 
field (BMF) it is possible to answer the question whether fragmentation is energetically 
favorable or not. In particular, a repulsive BEC in an asymmetric double- well potential can 
be in a stable two-fold fragmented state j3| , but the energy of this state is higher than the 
energy of the respective non- fragmented ground state of the condensate, i.e., it is an excited 
state of the BEC. This result is in agreement with general predictions derived by Nozieres 
P] that interaction prevents fragmentation in repulsive condensates. We shall demonstrate 
below, however, that in appropriate traps, repulsive BECs can exhibit fragmentation in the 
ground state. 

In the present paper we use the best mean-field approach to study the ground state of 
a system of N identical bosons with positive scattering length immersed into a double-trap 
external potential. A double-trap potential consists of an inner trap embedded in a wider 
outer trap. The paper is organized as follows. In Sec. II we briefly discuss the three-orbital 
best mean-field formalism. We also provide a very transparent example which illustrates 
that a finite number of particles and well-separated multiple potential wells are favorable 
conditions for fragmentation. In Sec. III. we propose specific shapes of trap potentials and 
discuss why fragmentation is expected. We demonstrate in Sec. IV that indeed in the ground 
state macroscopic occupation of three single-particle functions is energetically more favorable 
than accumulation of all particles in a single orbital. We also show that fragmentation is 
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a general phenomenon which, depending on the trap potentials, may take place for any 
number of particles. In Sec.V we demonstrate that the number of fragments, the shapes 
of the respective one-particle orbitals and their occupation numbers can be manipulated 
by the proper choice of the outer trap. Here, we also verify that fragmentation of the 
ground state is a critical phenomenon, because it occurs when the number of bosons exceeds 
some critical number at a fixed scattering length, or at some critical scattering length if the 
number of bosons is fixed. The interplay between critical parameters of the fragmentation 
and the geometry of the inner trap forms the content of Sec. VI. A discussion of the factors 
suppressing fragmentation in the ground state is given in Sec. VII. Next, we briefly address 
many-fold fragmentation in a multiple wells in Sec. VIII. Finally, Sec. IX summarizes our 
results and conclusions. 



II. MEAN-FIELD THEORIES 



We consider a system of N identical bosons interacting via a 5-function contact potential 
W(fi — fj) = XoS(fi — fj), where f\ is the position of the i-th boson and the nonlinear 
parameter Ao is related to the s-wave scattering length of the bosons j"!^ . 

The standart mean-field description of the interacting system is obtained by assum- 
ing the ground state wave function $ to be a product of identical spatial orbitals (p: 
r 2 , . . . , fjv) = (p(fi)(p(f2) • ■ ■ <p{tn)- The energy E = < >, defined as the 

expectation value of the H, reads 

E GP = N{J ^*h ¥ >df+^ J \^df}, (1) 

where A = Ao(A — 1) is the interaction parameter and h(r) = T + V(r) is the one-particle 
Hamiltonian consisting of the kinetic operator T and the external potential V(f). The mean- 
field equation for bosons residing in a sin gle orbital obtained by minimizing this energy is 




the well-known Gross-Pitaevskii equation \21\ . [22 

{ h(r) + A (A - l)Mr)| 2 } = Mgp (2) 

By definition, this equation cannot describe fragmentation because all bosons reside in a 
single orbital. The reduced one-particle density operator and the corresponding spatial 
density are given by PGp(f,r') = (p*(r')ip(r) and Pgp{t) = Iv 9 ^)! 2 respectively. 
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In order to describe the fragmentation on a mean-field level, a more general ansatz for 
wave-function must be used. In the present study three orbitals 0i, 2 and 3 with particle 
occupations ni,n 2 , n 3 and n x + n 2 + n 3 = iV are utilized. It will become evident below why 
three orbitals have to be used in the present investigation. Since the details of the general 
derivation have already been published elsewhere [^(J , we outline here only the main 
steps. With this ansatz the wavefunction now reads 

. . . , fjv) = <S0i(fi) • ■ ■ 0i(f ni )0 2 (r m +i) • • ■ M^M^m+i) ■ ■ ■ <fi 3 (r ni+n2+n3 ), (3) 
where S is the symmetrizing operator. The energy expression takes the form: 

g = n 1 / t u + A Wl(Wl 2 " 1) / \^dr + n 2 h 22 + X, n2{n2 - l) f \^dr+ 
n 3 h 33 + Ao ^ 3 ^ — — J |0 3 | 4 df + +2A nin 2 J |0i| 2 |0 2 | 2 df+ 

+2A n 1 n 3 ^ |0 1 | 2 |0 3 | 2 rff + +2A n 2 n 3 / |0 2 | 2 |0 3 | 2 df (4) 

By minimizing this energy with respect to the orbitals under the constraints that they are 
orthogonal and normalized, i.e., < <pi\(pj >= we get the following three coupled equations 
for the optimal orbitals: 

{ h(r) + A Q (m - l)|0i (f)| 2 + 2A o n 2 |0 2 (f)| 2 + 2A o n 3 |0 3 (r)| 2 } x (r) = 

= <j>i{r) + A*12 02 {r) + ^13 fo(f) 
{ h(f) + X (n 2 - l)|0 2 (f)| 2 + 2A o n 1 |0 1 (f)| 2 + 2A o n 3 |0 3 (r-)| 2 } 2 (f) = 

= A*21 01 (r) + ^22 02 (r) + /U 23 3 (f) 

{ h(r) + X (n 3 - l)|0 3 (f)| 2 + 2A o n 1 |0 1 (f)| 2 + 2A o n 2 |0 2 (f)| 2 } 3 (f) = 

= 1^31 M^> + ^32 M^) + ^3 Mr)- ( 5 ) 

Several key features of this approach should be mentioned. First, this mean- field in- 
cludes GP as a special case: when the occupation of two of the three orbitals vanish, 
i.e. n 2 = n 3 = 0, the system of equations (jSJ is reduced to the single equation (J2J) and 
the respective energy in EqE] coincides with the GP one (compare Eq0] and EqGJ). Sec- 
ond, by construction, this method can describe fragmentation for a finite number of frag- 
ments. Indeed, in this three-orbital case the reduced one-particle density operator can be 
written as Pbmf(Ti '"') — Yli 4>i { r ')4>i{r) n i an d the corresponding spatial density becomes 
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Pbmf(t) = ni\(f)i(r)\ 2 + n 2 |0 2 (^*)| 2 + ^3|03(^*)| 2 - Third, the occupation number n« of all frag- 
ments are variational parameters, minimizing the full energy. In order to find this optimal 
value of the energy, Eqs|5]are solved for all possible occupation numbers, and the respective 
energies, Eq0J are evaluated and compared. As we discussed before, the results obtained 
for the specific occupation numbers n 2 = n% = are identical to the standart GP ones. 
Therefore, within this computational scheme we automatically clarify the question on the 
favorability of fragmentation. Fourthly, A = Xq(N — 1) is the only parameter involved in the 
GP energy per particle expression (see Eqfl]). Therefore, the one-orbital mean-field cannot 
distinguish energetically between two system made of different numbers of bosons if they 
are characterized by the same A. In contrast to that, BMF treats these systems differently. 
For example, at the GP mean-field level two systems made of 11 bosons with Ao = 0.1 and 
of 1000001 bosons with Ao = 0.000001 are characterized by the same energy per particle. At 
the BMF level of description, if one of the systems is fragmented, the respective energies are 
different due to different fractional occupation numbers. This observation defines A, N as a 
very natural choice of the parameters to study fragmentation at the general mean-field level. 
If fragmentation takes place, then the number of particles becomes a relevant parameter. In 
the following, to compare systems made of a different number of bosons we adjust their Ao's 
in such a way that A is the same for each system. Then, for a given trap potential the GP 
ansatz gives the same energy per particle for all these systems. 

Large numbers of bosons and well-separated quantum levels of the single- well trap poten- 
tial provided a justification of the single-orbital mean-field description of BEG However, 
recent experiments on optical trapping of BEC have initiated interest in studies of multi-well 
systems and questioned the validity of the one-orbital mean-field description. The ground 
and lowest excited states of a multi-well trap potential can be almost degenerate and this 
opens a competition among the involved single-particle levels and the degree of their oc- 
cupations. Therefore, the one-orbital mean-field description may be insufficient. Another 
difficulty for the one-orbital mean-field arises when the number of wells is comparable to 
the number of particles. In this case the average occupation number of each well (so called 
filling factor per lattice site ) can be of the order of several atoms. It is worthwhile to 
demonstrate that for such a situation a many-orbital mean-field is the best mean-field, since 
it is energetically and physically more favorable than the one-orbital ansatz. 

Let us consider a system of 3 repulsive bosons trapped in three equivalent, infinitely 
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separated wells which we denote left (1) , central(c) and right (r). Without loss of generality, 
we can assume that the lowest total energy is obtained if each well contains only one particle, 
i.e., rii = n c = n r = 1. Since the wells are equivalent and infinitely separated, the wave- 
functions (orbitals) of each boson ipi, ip c , ip r are the same but localized at the different wells. 
This implies also zero overlap between each pair of orbitals. In other words, we have a 
system of three non-interacting bosons. The GP orbital now reads cp = (tpi + ip c + ip r ) / V3 
while the BMF orbitals are the three orbitals <pi = ipi, 02 — and 03 = ip r . A substitution 
of these functions into EqUJand EqHJ gives: 

E 1 - orh = h u + h cc + h rr + 3{3 ~ 2 1)Xo |(I|^| 4 + I|^ c | 4 + I|V; r | 4 )rff 
E^-orb = hu + h cc + h rr 

Inspection of these energies clearly shows that for repulsive interaction Ao > the three- 
orbital description is energetically more favorable than the one-orbital (GP) approach. The 
total energy of three equivalent non-interacting bosons is expected to be a sum of the energies 
of each particle. Therefore, only the BMF describes the physics correctly, while the GP 
energy contains an artificial term which can be considered as an interaction between the 
actually non-interacting subsystems. For this example of three infinitely separated wells, 
the energy provided by the many-orbital mean-field (BMF) is lower than the GP one for 
any finite number of particles. At the limit of very large N the energy difference between 
GP and BMF vanishes while the physics in the GP case still remains wrong. For this trivial 
case the three-orbital ansatz is evident. However, the question whether the number of the 
fragments is always equal to the number of wells, deserves a more detailed investigation and 
will be reported elsewhere j^j. 

We close this section with a general remark concerning BMF. For a given number of 
orbitals m , three in the present case, the BMF approach determines their optimal occupa- 
tion numbers minimizing the energy functional in EqHJ We call the result BMF(m ). As 
mentioned above, the calculation may provide the result that some of these orbitals are not 
occupied, i.e., their occupation number is equal to zero. This, of course, implies that the 
overall best mean-field is achieved with less than mo orbitals. Generally, we arrive at the 
overall best mean-field if inclusion of more orbitals does not improve the description. In 
the present study m = 3 and we have arguments that this choice leads to the overall best 
mean-field. 
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III. PROPOSED DOUBLE-TRAP POTENTIAL 



Our proposed ID double-trap potentials is shown in FigO Effectively such a trap may 
be obtained as a superposition of two potentials (inner and outer). We model the inner 
potential as: 



where A and B are parameters of the inner trap. As an outer trap embedding the inner one 
we used either an infinitely deep square potential well (infinite square well) with half-width 
equal C or a smooth power potential V ou t er (x) = (0.035a;) 10 . The infinite square well is 
obtained by placing the infinite walls at [— (C + 8) : +{C — 8)]. We introduced a small 
asymmetry parameter 5 = 0.017T to destroy the exact symmetry of the trap potential in 
order to get rid of the effects of symmetry. In the following, as a default outer trap we use 
a square well with walls at Cq = 9.5tt and A = 0.8, B = 0.1 as reference parameters of the 
inner trap. The corresponding kinetic energy reads T = — ^ -B^ implying that coordinate x 
is dimensionless and all energies and Ao are now in units of the frequency u. 

An inner tra p p otential of the form given in EqlHlhas been originally proposed to study 
BEC tunneling 0,13 • For the above reference parameters and one particle this potential 
has only a single bound state and a set of metastable states at positive energies (so called 
resonances). In the limit of non-interacting particles, all bosons will occupy this bound 
orbital, while for a non-zero repulsive interaction the bosons trapped inside the well may 
"flow out". In this case, a competition between the bound state localized inside the well 
and the continuum states outside the barriers has been predicted [25]. The experimental 
observation of the continuum outgoing waves is a delicate problem. By placing a secondary 
trap potential beyond the barriers the continuum outgoing wave functions are discretized 
and "transformed" to real functions which can be occupied by bosons and observed experi- 
mentally. We shall demonstrate later that the specific shape of the secondary trap potential 
is rather of minor importance, while its width is a major factor. The secondary trap poten- 
tial should have a width capable to accumulate particles. The simplest choice is to place 
two infinite walls at some distance from the origin. The fixed infinite walls can be replaced 
by some external embedding potentials with a smooth profile. From an experimental point 
of view this means that the fragmentation phenomenon may be observed in outer ordinary 
traps. 




(6) 
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IV. FRAGMENTATION 



The combined potential {Vi nner + V ou t er ) has three well-separated wells, see FigJTJ There- 
fore, if fragmentation takes place, bosons will be accumulated in each of these three wells. 
More precisely, the reduced one-particle density of the system of N identical bosons in this 
double-trap potential would have three macroscopic (with respect to N) eigenvalues. The 
condition that all three wells are well separated from each other implies that the respective 
eigenvectors (natural orbital) will be predominantly localized in each of these wells. 

Indeed, within the framework of BMF(3) the ground state of the system of N bosons 
becomes three-fold fragmented. In Fig|2] we present two sets of the orthonormal BMF 
orbitals and the respective densities corresponding to the system of N = 25 and N = 6000 
bosons and compare them with the GP results. The interaction strengths of these systems 
have been chosen to keep the quantity A = X (N — 1) = 1.3 fixed for both systems as 
explained above. BMF(3) predicts fragmentation of the ground state as it provides a lower 
energy than GP does. 

From FigJ21 one can see that the spatial densities (not to be confused with the reduced 
one-particle density) of GP, pep = \ f\ 2 , and of the three-orbital BMF, Pbmf = ^i|0i| 2 + 
^2 1 02 1 2 + ^3|03| 2 , are rather similar. At the same time the energy per particle provided by 
the BMF is lower than the respective GP one. Despite the fact that the spatial densities 
provided by GP and BMF are similar, there is a substantial physical difference between the 
methods. At the 1-orbital level of description the systems are unfragmented, while at the 
3-orbital one they are three-fold fragmented. 

By solving Eqs|S] for different occupation patterns we obtained the optimal fractional 
occupation numbers for the systems of N = 6000 and N = 25 bosons. For the sake of 
convenience we will from now on also use the term fractional occupation number for rii/N 
and express this quantity in % (^-100%). These are found to be similar for both systems: 
ni/N pa 68.58% for the orbital localized in the central well and n,2/N pa n^/N = 17.71% 
for the orbitals localized in the outer wells for N = 6000 and n\/N ~ 67.5% and U2/N pa 
n 3 /N = 16.25% for N = 25. The fact that the occupation numbers of the orbitals localized 
in the left and right wells are approximately the same is explained by the slight asymmetry 
of the double-trap potential. This observation allowed us to simplify the numerical search 
for the optimal values of the occupation numbers. Instead of searching for the minimum of 
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a functional E(nx,ri2, TI3) of two independent variables (723 = N — n\— 722), we can start the 
search using n-i = n.3 and then relax this condition. 

It is convenient to use the fractional occupation n\/N of the orbital localized in the 
central well as a characteristic parameter of fragmentation. If the fractional occupation 
n\/N — 100% then there is no fragmentation at all, while for any other values of n\/N the 
fraction of bosons accumulating in the outer wells is defined as 100%— n\/N . In the following 
we call the latter quantity the fragmented fraction in the outer wells or briefly the fragmented 
fraction. For the specific examples depicted in Fig|2 (A = 1.3), the fragmented fractions of 
the systems with N = 25 and N = 6000 bosons are 32.5% and 31.42%, respectively. 

This observation allows us to conclude that if the ground sates of the system made of a 
large number of bosons is fragmented, then any other system of bosons characterized by the 
same A and made of a smaller number of particles (of course, N > 1) is also fragmented. 
The opposite does not apply, however. For a given value of A there is a maximal number 
of bosons for which the ground state is fragmented. This number depends on the trap 
potentials used and can be manipulated by changing these potentials. For a discussion of 
this issue, see Sec. VII. 

In the following study we confine ourselves to the system of N = 25 particles, keeping 
in mind that for a larger number of particles the occupation numbers may differ within less 
than 5% as long as the condensate is fragmented (this has been verified numerically). 

V. MANIPULATING FRAGMENTATION BY VARYING THE OUTER TRAP 

By "manipulation of the fragmentation" we mean the possibility to choose the shape of 
the trap potential as well as the number of bosons (and possibly also their scattering length) 
in such a way that all fragments acquire the desired occupation numbers. 

As we briefly mentioned above, the bosons trapped by the inner potential alone in the 
absence of the outer trap occupy the bound state (localized in the central well) as long 
as A < A cr . Any change in the number of particles or in the scattering length such that 
A = X (N—1) becomes larger than A cr , imidiatelly initiates tunneling - the flow of bosons out 
of the central well [3] . If infinite walls are placed beyond the barriers, the system becomes 
closed and bosons are collected in the outer wells. If all three wells are macroscopically 
occupied the system is fragmented. 
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Several questions arise in the presence of the outer trap. Does the fragmentation phe- 
nomenon exist for any A or is it characterized by some critical parameters, similarly to 
tunneling in the open system? The second question may be formulated as follows: Do the 
fractional occupation numbers depend upon the positions ±C of the walls of the outer trap 
and on the particular shape of this trap? 

Figure H3 shows the fractional occupation of the orbital localized in the central well as 
a function of the positions ±C of the outer walls (see Figure for several values of A. 
From this figure it is clear that the fragmentation starts to take place when A exceeds some 
threshold (for the example of 25 bosons and C = lln X cr = 0.8249). It is interesting to 
notice that the exact value of this threshold for a system with a finite number of bosons is 
slightly smaller than that for the open system, i.e. C — > oo. As N grows, the critical value 
of A obtained for the closed system approaches the numerical result A cr = 0.8279 for the 
open one where tunneling through the barriers begins 2^. 

Further increasing the bosons interaction strength A beyond A cr (at least up to A = 3.0) 
leads to a more pronounced fragmentation of the ground state. Here, we have to mention 
that there is another limit where fragmentation must disappear, namely when A becomes 
so large that the chemical potential is larger than the barrier heights of the inner trap and 
particles can flow freely into the outer trap. 

These observations reveal that fragmentation of the ground state is a critical phenomenon 
initiated when the number of bosons exceeds some critical number, for a fixed scattering 
length, or at some critical scattering length if the number of bosons is fixed. 

Figure El also illustrates that fragmentation can be observed if the width of the outer 
trap is quite large or more specifically, if the infinite walls are placed at comperativelly large 
distances from the barriers. For example, for the chosen shape of the inner trap (Aq, Bq) 
with the barriers at x^ ~ ±3.4, and A = 1.3, about 30% of bosons are shared by outer wells 
if the walls are at C ~ 25. By pushing the walls toward the barriers, the fragmentation 
gradually decreases: for C ~ 16 the fragmentation is 20%, and for walls at C < 10 the 
fragmentation totally disappears - all bosons are accumulated in the inner trap. Therefore, 
by pushing the walls toward the barriers the fragmentation can be reduced or even totally 
suppressed. These results show that fragmentation in the double-trap potential may be 
manipulated by the outer trap: fragmentation can be suppressed by squeezing the outer trap 
and enhanced by its expansion. 
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On the other hand, when the walls are moving outward from the barriers, the fragmented 
fraction in the outer wells becomes larger and converges to some constant value which 
depends on A, of course. This value can be extracted from the results of Fig |H] and also from 
the results obtained for the open system; an open system may be thought of as a closed one 
with walls placed at infinity. From this we conclude that for a given value of A, the maximal 
fraction of the bosons in the outer wells is defined by the inner-trap potential only. 

The relevant factor for the fragmentation is the width of the outer trap and we may 
suppose that the specific shape of the outer trap is of lesser importance. To support this 
expectation, we study the fragmentation for A = 1.3 in the systems of N = 3000 and 
N = 25 bosons trapped in the double-trap, with the same inner potential as before and 
the smooth power (0.035x) 10 outer potential (red solid line in Fig. [TJ. As seen in Fig0J 
the three-fold fragmentation of the ground state is again favorable energetically for these 
systems. Moreover, the optimal occupations of the inner orbital for iV = 3000 and N = 25 
are rii/N 72.2% and n\/N ~ 71% respectively, and hence similar to those discussed 
above for the infinite walls case. In Fig. 0] we plotted the orbitals and the spatial densities 
determined for the smooth power outer trap. By comparing these orbitals with those shown 
in FigOfor the infinite square outer trap, we conclude that the shape of the orbital localized 
in the inner trap does not depend upon the specific shape of the outer potential. The profiles 
of the orbitals localized in the outer wells exhibit differences which may be experimentally 
observed. In the case of the smooth power outer trap the density profiles of the outer orbitals 
are gaussian-like, whereas those for the infinite square outer trap are of a sinusoidal type. 

VI. MANIPULATING FRAGMENTATION BY VARYING THE INNER TRAP 

Here, we investigate how fragmentation of BEC depends upon the shape of the inner 
potential. This potential permits two degrees of manipulation by varying the depth of the 
inner well and by varying the height of the barriers. The parameter A of the inner potential 
(see Eqinjand Fig|SjA.) is directly related to the depth of the inner well. This depth grows 
as A is increased. In Fig. 03 we illustrate the potential dependence upon B, a parameter 
which defines the height of the barriers, their widths, and the positions of their maxima. 
By decreasing the value of B, the height of the barriers and their widths grow while the 
positions of their maxima are shifted outward. 
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In Fig El we plot the fractional occupation in the inner trap as a function of A for several 
values of A. It is seen that fragmentation decreases monotonically with A. Qualitatively 
speaking, this finding implies that by increasing the depth of the inner potential more "room" 
becomes available for bosons in the inner well. Conversely, by decreasing A the capacity 
of the inner well becomes smaller and more particles "flow" out into the outer wells. This 
picture serves also as a verification and extension of the conclusions drawn above. We 
have established that at fixed trap geometries fragmentation takes place at A ~ A cr and 
becomes more pronounced as A is increased. We may now add that there is a critical depth 
of the inner potential initiating fragmentation. For example, for A = 0.8 and A = 0.8 
(B = 0.01, C = 9.57r) fragmentation does not exist, while for A = 0.72 the system becomes 
20% fragmented. 

Finally, we investigated the dependence of the fragmentation on the height of the barriers 
(5-parameter). In Fig. Q we plot the fractional occupation of the inner well as a function 
of B (A = 0.8, Co = 9.57r) for several A > A cr . From this figure it is clear that as long as 
B is not too small or too large the fragmentation is not particularly sensitive to variations 
of B. Decreasing B from 0.16 to 0.02 corresponds to a substantial change of the barriers 
heights from 0.8 to 5.8 units. At the same time the respective fractional occupation of the 
bosons in the central well varies by several percents only for any fixed value of A. A further 
decrease of B causes, however, the disappearance of the fragmentation. Such a behavior is 
to be expected in this case, since a very small value of B corresponds to very broad barriers 
(see FigEJ) diminishing thereby the size of the outer wells and hence their capacity to hold 
bosons at favorable energy cost (see. Sec.V). 

On the other hand, at large values of B the height of the barriers becomes very small 
and the fragmentation of the ground state is expected to be unfavorable energetically. This 
issue is further discussed in the subsequent section. 

It is very important to note that the presence of barriers is an essential factor for the 
ground state fragmentation of the repulsive BEC. The specific shape of the barriers is of 
lesser relevance and its impact on fragmentation can be largely compensated by moving the 
position of the walls of the outer trap. We expect, however, that for time-dependent studies 
this situation will be changed drastically, because the tunneling time (i.e. the time which is 
needed to tunnel through the barriers) is determined by the height and width of the barriers. 
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VII. MANIPULATING THE MAXIMAL NUMBER OF BOSONS IN THE FRAG- 
MENTED GROUND STATE 



There are several factors limiting the number of bosons in a fragmented ground state. 
The most transparent factor is the height of the barriers of the inner trap. Increasing the 
particle number N obviously enlarges A = Xq(N — 1) for a given scattering length. We have 
already pointed out that when A becomes so large that the chemical potential is larger than 
the barrier heights, particles can flow freely out of the inner trap and the fragmentation 
disappears. Consequently, there is a maximal number N max of bosons in a fragmented 
ground state and this number depends on the double-trap potential and on the scattering 
length. 

For the open system we have seen that- at a fixed value of A- once fragmentation in 
the ground state takes place, i.e., A > A cr , this fragmentation persists for any number of 
bosons (N larger than 1, of course). This finding does not hold for closed systems. For these 
systems there is again an N max even if A is kept fixed at a value where the chemical potential 
is smaller than the barriers heights. As seen in the preceding sections, there is an enormous 
range of double-trap potentials giving rise to fragmented ground states as long as A > A cr . 
The degree of fragmentation can be widely manipulated by varying the parameters of the 
trap potentials. Two major questions arise now: what is the origin of N max and can this 
value be manipulated by varying the parameters of the trap potentials? The answer to the 
latter question is positive and we shall present numerical calculations below. 

To facilitate the origin of the existence of N max we remind that A = X (N — 1) where A is 
the interaction strength of two interacting bosons. Increasing the number N of bosons while 
keeping the value of A fixed obviously implies a weakening of the interaction strength. As 
N approaches infinity, the interaction strength Ao approaches zero and we may expect that 
the ground state takes on the appearance of a non-interacting system, i.e., that of the GP 
ansatz. Indeed, our numerical calculations show that while for N < N max the ground state 
is fragmented, fragmentation still persists even for N > N max , but the state in question is 
an excited state of the condensate. This excited state is very low lying for all N > N max ; 
its energy is extremely close to that of the now unfragmented ground state (in the present 
examples the two energies per particle differ just at the 7 th digit!). The existence of a very 
low-lying fragmented excited state is of great interest by itself and should play a role in 
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particular in time- dependent experiments. 

In Fig|Hlwe show N max as a function of the outer trap parameter C for several values of 
A keeping the parameters of the inner trap at their reference values. Also shown are the 
variations of N max with the inner trap parameter B. In the range of parameters studied, 
the maximal number of bosons N max grows linearly with the size of the outer trap and 
exponentially with the height and width of the barriers. 

For a better understanding of the above findings we draw attention to FigO To maintain 
orthogonality, the three orbitals of the BMF live in different regions of space. As long as 
the size of the outer traps is large enough, this is favored by the presence of broad and high 
barriers between the inner and outer wells. Since the condensate is repulsive, a large outer 
trap enables the wavefunction in this region of space to delocalize without the necessity to 
penetrate the inner well. 

VIII. MANIPULATION OF FRAGMENTATION BY SEVERAL INNER TRAPS 

Three-fold fragmentation of the ground state is found to take place in suitable three-well 
potentials. The results obtained for these potentials can naturally be extended to an array of 
multiple wells. This extension is particularly relevant due to the recent experimental efforts 
to investigate one-dimensional optical lattices (12j. For example, a multi-well potential may 
be formed if the inner- well trap given by EqGDis translated several times and then embraced 
by an outer smooth power or infinite square trap. In these cases we expect many-fold 
fragmentation and if the inner traps are well separated from each other and from the outer 
walls, then the critical value of A scales according to the number of inner potentials. In FigEl 
we plotted an example of such a trap. The critical value of A for this system is approximately 
five times larger than the respective value obtained above for the single inner potential. This 
value as well as the extend of the fragmentation can be sensitively manipulated by varying 
the trap parameters. Of course, a many-orbital BMF must be used in order to correctly 
describe the situation. For the multi-well trap shown in Fig|U]we expect that 11 orbitals 
should be used and that should be possible to do in the future. 
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IX. CONCLUSIONS 



In this article we have investigated the fragmentation phenomenon in the ground state 
of a repulsive condensate immersed into the double-trap potential. We demonstrate that 
fragmentation can be successfully characterized by the best mean-field approach. To be 
able to correctly describe m-fold fragmentation, m orbital are needed and are available in 
the m-orbital mean-field BMF(m) method. In this method the occupation number of each 
fragment as well as the optimal shape of the respective orbitals are determined variationally 
by minimizing the total energy functional. If more orbitals are included than needed, the 
occupation of the superfluous orbitals becomes zero by minimizing the total energy. 

The double-trap potential studied here has three wells, obtained as the superposition of 
an inner trap exhibiting one well and two barriers and an outer trap embedding the inner 
one. For many choices of the potentials, the macroscopic occupation of the three orbitals 
may become energetically more favorable than accumulating all the particles in a single 
orbital. The fragmentation of the ground state is found to occur when the number of bosons 
exceeds some critical value which depends on the scattering length and on the shape of the 
inner trap potential. For the example studied we found that if fragmentation is observed 
for a large number of bosons, then it exists also for any smaller number of bosons (of course 
N > 1) when A is kept fixed. When A is kept fixed, there exists, however, a maximal 
number of bosons for which the ground state is fragmented. This number can be strongly 
manipulated by varying the double-trap potential. 

We have demonstrated that the geometry of the inner potential determines the values 
of the critical parameters. Moreover, for any given number of bosons this potential also 
determines the maximally possible fragmented fraction of bosons which is localized in the 
outer wells. The actual fragmented fraction of bosons may be effectively manipulated by the 
proper choice of the outer trap. The interplay between the inner and outer trap potentials 
provides a sensitive tool to manipulate fragmentation of repulsive condensates. Varying the 
number of bosons in the condensate and the scattering length are also instrumental in this 
respect. 

The results obtained for three-well potentials can naturally be extended to an array of 
multiple wells. 
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FIG. 1: Proposed double-trap potential made of an inner trap of the form given in Eq. Eland outer 
trap. As outer traps we use either infinite walls at ±C (black) or a smooth power potential (red). 
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FIG. 2: The orbitals and densities for the double-trap potential with infinite walls (see black curve 
in Fig^). The orbitals cj>i(x), 02 0*0 and <t>$(x) and the respective density per particle pbmf = 
(tt-i |0i | 2 + ^2 1 02 1 2 + n3\(j)3\ 2 )/N of the three-fold fragmented ground state (for convenience the 
(ni/N) 1 / 2 (f>i(x) are shown) are depicted in comparison with the GP orbital <p and its density 
Pgp = \ip\ 2 for N = 25 and N = 6000 and A = 1.3. The energy per particle is indicated. For 
convenience, the base-line of the orbitals 4>2(x) and 03 (x) has been moved upwards artificially from 
zero. 
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FIG. 3: Manipulating the fragmentation by varying the outer trap. Shown is the fractional 
occupation of the orbital localized in the inner well as a function of C (half-width of the outer 
trap). All other parameters are kept at their reference values Aq and Bq. 
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FIG. 4: The orbitals and densities for the double-trap potential with a smooth power outer trap 
(see red curve in Fig^). The orbitals 4>i(x), 4>2(%) and ^{x) and the respective density per particle 
Pbmf = (ni\4>i\ 2 + ^Ifol 2 + n 3\4>3\ 2 ) / N of the three- fold fragmented ground state (for convenience 
the (rii/N) 1 / 2 4>i(x) are shown) are depicted in comparison with the GP orbital ip and its density 
Pgp = \<fi\ 2 for N = 25 and N = 3000 and A = 1.3. The energy per particle is indicated. For 
convenience, the base-line of the orbitals 4>2(x) and 4>3(x) has been moved upwards artificially from 
zero. 
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FIG. 6: Manipulating the fragmentation by varying the inner trap. Shown is the fractional 

occupation of the orbital localized in the inner well as a function of A. All other parameters are 

kept at their reference values Bq and Cq. 
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FIG. 7: Manipulating the fragmentation by varying the inner trap. Shown is the fractional 

occupation of the orbital localized in the inner well as a function of B. All other parameters are 

kept at their reference values Aq and Cq. 
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FIG. 8: Manipulating the maximal number of bosons N max in the fragmented ground state by 
varying the double-trap potential. Left figure: Dependence on the parameter C. All other param- 
eters are kept at their reference values Aq and Bq. Right figure: Dependence on the parameter B 
(note the logarithmic scale). All other parameters are kept at their reference values Aq and Cq. 
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FIG. 9: Proposed multi-well double-trap potential. 
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